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APPENDIX: FIGURES AND TABLES 
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The Lightning Detection and Ranging (LDAR) System developed over a number of years by 
Carl Lennon and colleagues at the Kennedy Space Center is based on the times of arrival at 
seven locations of electromagnetic radiation in the 66MH frequency range emitted by electric 
charge movement in thunder storms. The differences in the times of arrivals are converted into 
differences in distances from the point of origin (x, y, z ) of the radiation to the various receiving 
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we see that 



= 299792.458 km/sec 
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We note here that equations (2) for i, j* 0 do not admit this simplification and remain ncsalinear 
in x, y and r . The pairs of three equations i = 1, 3, 5 and i = 2, 4, 6 are solved for x, y aoc. ~ with 
z determined from 
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On the other hand if (12) is not satisfied, then each of the 20 combinations of three in lecendent 
equations from the six equations in (9) are solved and each solution is compared with ihe others 
to find at least eight solutions which satisfy the tolerance limits in (12) with R in (14) :u_>_:iilated 
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Given any three equations /, j, k, i* j * k, we can write 




Note that we have displayed the dependence of x, y, and r on i, j, k via subscripts. For perfect 
data, the x ijt , y ijk , and r ijk are all equal respectively. Suppose that u, is replaced by u t + £, and that 
v, is replaced by v f - u,£, - (ef /2) in (16) above. Then, the system (18) can be written as 


where the solution s iit denoted in (20) above has its dependence on the errors £,.,£., and £, 
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Subtracting (18) from (22), we see that 



We note two things. First, for fixed measurement errors, the overall solution error 
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For derivations and discussions of vector norms, matrix norms, condition number, and relative 
error, see [18] or any other reasonable text on numerical analysis or numerical linear algebra. For 
the vector norm defined as the maximum absolute value of the components of the vector, the 
corresponding matrix norm is the maximum of the row sums of the absolute values of the entries 
in each row; i.e., for the matrix A = ) 
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While GDOP is a problem for (x, y, r), it is a major problem for z. Consider the expansion 
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From spherical coordinates, we recall that z = rsin <f> where <p is the angle of elevation out of the 
at - plane. By differentiating the equations (22) with respect to £ i ,e j ,e k , we observe (Not 

Proved-although a messy induction argument looks feasible) that 
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electronics is subject to various interference sources that corrupt various portions of the data. 
Given the sensitivity of the system to errors, it is feasible that an error in the data could occur 
which would fake an acceptable solution to the system. It would be interesting to investigate the 
various possible types of errors in the system which could produce such fakes and whether or not 
the error would be an occasional fake or a frequent one. We shall comment on this again. 
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±.5 kilometer one could always find an error combination which would produce a fake solution, 
i.e., the solution corresponding to the least /i would be several kilometers away from the 
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all of the slack variables from the solution of the programming problem. Hence, we had to resort 
to a brute force approach of examining all solutions of combinations of three equations selected 
from the six equations in (16). For each such solution (x, y, r), we defined the residuals 


Then we ordered the residuals 




yields the solution (*, y, r ) corresponding to the best fit of five of the six equations in (16) 
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and the Newton method failed to converge. Gradient procedures to find a good starting set of u ‘s 
were tried and discarded as were simple search methods. Occasionally, the Newton method 
would converge to a fake solution especially when the linear procedure produced one. 
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E. Statistical Estimations. Lots of solutions can be generated when the data u ; ,i = l,—,6, 
contain error. For example, three equations from six yields 20 different solutions, four equations 
from six yields via least squares 15 more different solutions, five from 6 yields 6 more and least 
squares for all six yields another for a total of 42 different solutions which can be seen to 
coalesce into the correct (x, y, r) as the error in the data tends to zero. There are many ways to 
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Various data fit measures for each solution (jc t , y k , r k ) can be calculated. For example, 




solutions with the highest number of tally marks defined the cluster. The approximate solutions 
in the cluster could be viewed as reasonable approximations in a broad sense. However, as the 
data error increased fake solutions were observed. 

To filter out fake solutions some additional measurements could be very helpful. 
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. Conclusions and Open Questions 


00 ^ 
a O 


i 

4-» 

8 

I 

i 

I 

X) 

O 

O 

oo 

5 

• t -4 

o 

£ 


T3 

O 

•s 

<D 

a ® 

i ^ 

i < 

o 

U 


•o 

a 

p 

<2 

* 

c/5 

G 

0 
*G 
G 

1 

I 

00 

a 

a 

• »-4 

a 

• *H 

13 

U-. 

<2 

T3 

1 

<0 

S 

4 h 

0 

s 

Q< 

1 

O 

<2 

o 

2 


PQ 


G 


G 

O 

■a 


o 

S <§ 
o 00 
o o 

0 4J 

.2 2 
2 o 
2 2 

S- § 

<y .£ 

• S w> 

1 | 

I 8 

op ,8 

+3 £ 

g .a> 

a 3 


<o 


<y 

.52 

a 

2 


&I -a 

05 P 

£ 

o 


co O 

8 a 

.2 «3 

£ § 
00 *S 
S * 

G 


5 G 
O G 


43 

co 


a 

r— I 

E 

a 


<D 


CO 

2 

f 

a, 


a 

G 

<D 

_ > 
CD CD 

S SP 

g .g 

a> 

a 

o 

CO 


e s s 


u 


43 

00 


O 

*43 


<D 

a 

a 

a 

* 

4-* 

cx 

a 

1 

o 

2 

G 

O 

•G 

*»H 

G 

4h 

0) 

T) 


§ 

o 

•O 

<0 

<D 

G 

T3 

G 

O 

£ 


8 

m 




00 —4* 

<D T3 

*S 2 

<D 

T3 g 
O , 
G V* 
oo 

• • <U (D 

8 6 m 


G 

T3 


G 

a 3 a 
§ E* * 

g 

4> ^ 

04 

O 


Si 

3 

a 

<D 

a 


G 

2 

G 

£ 

co 

<0 

G 

(D 


G 

G 


PQ 


CO ^ 

1 2 

> CD 

2 <p 

G O 

• rH G 

CO S 

00 c 
O G 

<0 2 
2 2 
a ’a 

«» o 

a 2 
5 8 
£ w 

u 


<D 


8 

-3 

CO 8 

2 

c/5 2 

O CO 

G co 
w CD 

2 8 


00 

G 

• tH 

X) 

G 

O 

cx 

co 

8 

o 

o 

co 

<D 

CD 


G 

O 

*G 

G 

rH 

o 

CO 

<L> 

G 


G T 3 

6 *> 


>> 

co O 

2 "ex 

§ a 

G <D 
> 

<D 

G 


<U 

5 


1-4 

<D 

G 

6 

G 

G 

G 

O 

\G 

• rH 

T3 

G 

O 

CD 

<D 


co 

T3 

<1D 

<U 

G 


VO 

VI 

V 


G 
O 

=2 

v a 

VI & 


2 

a 

•a 

0 

•s 

1 

G 

a 

i 

O 

< 4-1 

O 

§) 

8 

Q 

ui 


o 



Cianos, N., Oetzel, G.N., and Pierce, E.T., A technique for accurately locating lightning at 
close ranges, J. Appl. Meteor., 11, 1120-1127, 1972. 


00 


d) 


2 

.s? 

* Kt 

s 

-G 

00 

•G 

I 

‘S 

ik 

<*> 

O 


Pd 

N 
O 
♦ *-< 

£ 

d> 

■9 

o 

00 


X 


•a 

•s 

o 

fcd 


c3 

o 

cl 

a- 


§ 

a t 

1 1 

u I 

r 

2 S 
* £ 

2 

§ 


I 


OX) 

.9 

s 

G 

OX) 


co 

<0 

*3 

3 

CO 

Uh 

d> 

T3 

G 

X 

c 

o 

• H 
+-» 
o 

d) 


PL, 

X 

> 

u 

od 


T) 

g 

G 

Q 


g 

<u 

<u 

0 

a 

s 


oo 

r- 

On 


uo 

On 

r- 

i 

vo 

oo 

r- 


k. 

o 


oo g 
g o 

5b 'S 


§ 


U. < 

G 

o 

d) 

a 

9 


§ 

a 
o 

g S 
JJ I 


<D 

Q 


o 

U 


a» 5 

•a .a 

1 ^ 
W) J> 

3 5 


00 


,<D c rr? 


U 

G 

O 

00 


X 


g 

o 

G 

G 


d> 

*3 

s 


co 

00 

.9 

T 3 

d) 

d> 

O 

o 


d § 
u 9 


CO 

G 

I 

a 

Q 

oo 

G 

X 

8 

S 

d) 

o 

o 

00 

73 

o 

•rH 

OX) 

O 


'O CL, 
§ co 

• rH 

r * "5 

- f 

•a 3 

<D 
O 


2 

o 

2 

<u 


*2 

£ 

a 

a 

CO 

>% 

00 

Pd 

< 

Q 

d 


cd wo 
rt O' 

SI © 

£ cs 
oo >o 


I ^ 

cd 3 
(D a 

£ I 


I 


G 

O 

$ 


OX) 

G 


G 

OX) 


X) 

G 

X 

G 

O 

X 

8 


PL 

as 

> 

Q 

£ 

G 
d) 

9 

d) 

0 g 

a ? 


G 

G 

G 


CO 

r- 

On 


O 

,, d- 


(J ^ 

j rH 

Pd . 

I § 

s 5 


O 

<N 

i 

ON 

ON 


d 

Co 

0 

1 


1 


Pi 

U 

ts 

o 

CL 

U 

Pd 


a 

X 

G 

O 

U 

< 

00 

< 


£ 


& 


ON 

a 2 


■S 




& 

2 

a 

o 

c 

H 

4 

ad 

<3 

3 

8 

CL 


On 


m 

vo 

H" 

IT) 


. 

3L On 
& & 


wo 

0 

1 -H 

Tt- 

r- 

i 

g 

a 

a 

T3 

1 

O 

S 

<u 

s 

‘I 

X! 

g 

< 

< 

z 

n ‘G 

>— I 
U 


I 

II 

M 

•r*» 

ftd 

R 

O 

’G 

S 

so 

Q 

do 

,R 
• *«■» 

s 

-R 

do 

•rH 


G 

O 

G 

G 


T3 

G 

G 


d> 

s 

d) 

o 

p^ 


d> 

o 

R 


d> 

CL 


< *8 
X 


R 

O 

*G 

.CL 

o 

^3 


oB 

G 

• rH 

3 

x< 

OX) 


CO 

2 

3 

o 

♦ rH 

CL 

O 

3 

G 

1-4 

Ph 

X 

> 

OX) 

G 

♦ rH 

G 

♦ rH 

G 
•*— > 

X) 

O 


(J) T— H 

a ^ 

>> — 
CO r ~‘ l 

O ON 
G OO 

o ^t 

x> "7* 

Vh 1 

S oo 
cl £: 

G ^ 

•N 

< s 


w 

d 




O 


§• 
PU O 


o 

o 


oo 

On 




Tj* 


to 

to 

§- 

tu 

o 


CO 

d) 

G 

CO 

G 

X 

-o 

G 

3 

73 


CO 

2 

a 

o 

• rH 

CL 

O 

3 

G 

Vh 

Ph 

X 

> 

pi 

d 


o 

o 


(N 


CO 


wo 


VO 


00 


On 



12. Proctor, D.E., Lightning and precipitation in a small multicellular thunderstorm, J. Geopkys. 



17. Twomey, S.A., Introduction to the Mathematics of Inversion in Remote Sensing and Indirect 
Measurements, Elsevier, New York, 1977. 


18. Young, D.M., and Gregory, R.T., A survey of Numerical Mathematics, Volumes I & II, 
Addison Wesley Publishing Company, Reading Massachusetts, 1972. 


» 


CM 





FIGURE 1: LOCATIONS OF LDAR SITES 




TABLE#1 : L-00 CONDITION NUMBER FOR THE MATRIX OF (1 9) 

WITH THETA= 45 DEGREES AND R=50 KM T 
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FIGURE 3: ANGLE FILTERED SOLUTIONS FROM ADMISSIBLE 
SOLUTIONS WITH MAXIMUM RANDOM ERROR 750 METERS AND 
SOLUTION AT X=15 KM, Y= 15 KM AND Z= 20 KM 
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FIGURE 5: ANGLE FILTERED SOLUTIONS FROM ADMISSIBLE 
SOLUTIONS WITH MAXIMUM RANDOM ERROR 750 METERS AND 
SOLUTION AT X=15 KM, Y=15 KM, AND Z= 20 KM 
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